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We construct a mass matrix for the four neutrino flavors, 
three active and one sterile, needed to fit oscillations in all 
three neutrino experiments: atmospheric, solar, and LSND, 
simultaneously. It organizes the neutrinos into two doublets 
whose central values are about 1 eV apart, and whose split- 
tings are of the order of 10 -3 eV. Atmospheric neutrino os- 
cillations are described as maximal mixing within the upper 
doublet, and solar as the same within the lower doublet. Then 
LSND is a weak transition from one doublet to the other. We 
comment on the Majorana versus Dirac nature of the active 
neutrinos and show that our mass matrix can be derived from 
an S2 x S2 permutation symmetry plus an equal splitting rule. 

PACS numbers: 14.60. Pq, 12.15.Ff 



Neutrinos produced by the interaction of cosmic rays 
with the Earth's upper atmosphere provide the strongest 
evidence for neutrino oscillations |IJ, with — > v T as 
the favored flavor transition If the additional evi- 
dence from solar ||] and LSND Q experiments is also 
confirmed, then it will be necessary to introduce a fourth 
light neutrino, a so-called "sterile neutrino" v s in addi- 
tion to the standard electron-, muon-, and tau- neutrinos 
to account for all the data ||. The question then arises 
as to the mass spectrum and mixing scheme for these 
four particles. 

In a two-flavor oscillation scenario, the atmospheric 
data suggest maximal mixing with mass difference 
Am 2 « 3 x 1CT 3 eV 2 [§. Of the three types of solu- 
tion for the solar neutrino data, there are two, namely 
the large angle MSW (LMSW) and the "just-so" in vacuo 
ones, which require close to maximal mixing 0; while the 
third, small angle MSW (SMSW) requires small mixing 
||. In all three cases, the mass difference Am 2 is much 
smaller than in the atmospheric case. By contrast, the 
LSND data require small mixing, but with a relatively 
large Am 2 as compared with the atmospheric case [[| . 

Recently, Bilenky, Giunti, Grimus, and Schwetz j^] 
have shown that the only way to account for these data 
in a four neutrino framework is to require a mass spec- 
trum consisting of two doublets |fTofl , with the splitting 
within each doublet being much smaller than the sepa- 
ration between them. Here we wish to propose a specific 
realization and mass matrix in which the members of the 



upper doublet are identified as maximal superpositions 
of and v T , and the members of the lower doublet are 
maximal superpositions of v e and v s . Atmospheric neu- 
trino data can then be described as maximal oscillations 
between the levels of the upper doublet, and solar neu- 
trino data as maximal oscillations between the levels of 
the lower doublet. LSND is then a weak transition from 
one doublet to the other. 

In adopting this point of view, we recognize that there 
are some problems with the current data. The valid- 
ity, or otherwise, of sterile neutrinos will be tested at 
SNO. [[llj The main impetus for reviving the "just-so" 
solutions comes from the anomalous points at the high 
end of the solar electron recoil spectrum observed at Su- 
perKamiokandc. [Q As we have pointed out in another 
paper, ||l2|| a crucial test for this will be the measurement 
of the 7 Be neutrinos. 

Our approach to the development of a mass matrix for 
a two-doublet model can be illustrated with the following 
two-dimensional model: 



m s m k \/tp a 



(1) 



in which the matrix M2 is a linear combination of the 
unit (2 x 2) matrix / and the Pauli matrix o~ x : 



M 2 = m s I + m k a x 



(2) 



It has eigenvalues [m s ± rrik) and eigenstates which are 
maximal mixtures of the basis states: 



(3) 



and thus it will lead to maximal mixing between neu- 
trinos v a and Vb- For future reference, we note that the 
matrix M2 is symmetric under the permutation group S2 
of the two members of the doublet, and that the eigen- 
vectors ip± are respectively even and odd representations 
ofS 2 . 

Now suppose we rotate M 2 through a small angle 
(—269) about the y-axis: 



exp (+ia y S6)M2 e^q>(—ia y S8) 
— m s I + vnkO~ x cos 259 + m k = o z sin 259 
rn s + rrik sin 259 — m k cos 259 \ 
rrik cos 259 m s — m k sin 259 I 



(4) 
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It has the same eigenvalues as the original matrix, but 
its eigenstates are also rotated through the small angle 
(-266): 



ip + cos 59 + tp_ sin 59 
ip + sin 59 — V- cos 



(5) 



and so it leads to small mixing oscillations between ip + 
and ip- . 

Guided by this analysis, we propose a four-flavor mass 
matrix which we construct by replacing m s and m k in 
the rotated form of M 2 by (2 x 2) matrices: 



m s -> M , M = 
m k —>K, K = 



m s m d 
m d m s 

m k 
m k 

Our model then takes the form 
*M 4 * = 

(V ^,(M + Ksm256 Kcos269 \{* a \ 
K*a *b)y Kcos 269 M-Ksm.256) J ' 

where ^ a and arc now two-dimensional column vec- 
tors: 



(6) 
(7) 

(8) 



_ "002 ^62 

Next we rotate M4 and 'J into the forms: 
,'M K 

Ma —>Ma = 



K M I ' 



(9) 



(10) 



$ fc y ^smW cos <56> y ^ * 6 y v 7 

For future reference, we note that M 4 is symmetric under 
the permutation group £2 which interchanges the two 
doublet pairs. 

Now we rotate M 4 and $ into: 



M 4 



M + K 
M-K 



y/2~\-<S>a + $6 



(12) 
(13) 



We now have to diagonalize the (2 x 2) matrices (M± 
if), where 



to s ± TOfc m d 

m d m s ±m k J ' 



(14) 



which have eigenstates ($ a + $b)/\/2 and 
(-$ a + $b)/\/2 respectively. The eigenvalues of (M+K) 
are: 



and those of (M — K) are: 

Mg =m s -m k ±m d 



(16) 



Thus we have two doublets whose mean masses are sep- 
arated by 2m k , and whose splittings are both given by 
2m d . The upper and lower components of ($ a + <!>{,)/ -s/2 
are maximally mixed, as are those of (-$ a + $ 6 )/-y/2. Fi- 
nally, the eigenstates of (M + K) are weakly mixed with 
those of (M — K) via the relation between $ and in 
Eq. 11 above. 

We identify (M + K) and its eigenstates with the at- 
mospheric neutrino oscillations between and v T , and 
so the squared mass difference may be written 

A A = (m s + m k + m d ) 2 - (m s + m k - m d ) 2 

= 4(m s + m k )m d . (17) 

Similarly, we identify (M — K) and its eigenstates with 
solar neutrino oscillations between v e and v s , and so 

A s = (m s -m k + m d ) 2 - (m s - m k - m d ) 2 

= 4(m s - m k )m d . (18) 

For reasons which will become apparent below, we write 



m fe = too - e 



and so 



e 

too 



As 
A A 



(19) 
(20) 



(21) 



Since A A is much greater than As, as discussed below, 
we conclude that e is much smaller than m , and that 
to s is only marginally greater than m k : 



m s 
m k 



(l + 2e) ~ 1 + 2 



As 
A A 



(22) 



For LSND, we assume that the — > v e oscillation is 
dominated by the transition from the lower eigenvalue of 
(M + K) to the upper eigenvalue of (M — K): 



A L = (m s + m k - m d ) 2 
= 4(m k - m d )m s , 



(m s - m k + m d ) 2 



and so 



8m k m s 



2 A, 



- £ 2 ) 
A A 



(23) 



(24) 



Since A^ is much bigger than either A^ or As, it follows 
that: 



2 to 



'A, 



1 + 



Aa 
4A, 



(25) 



M£ = m s + m k ±m d , 



(15) We then find that m d is much smaller than to : 



2 



2m d 



A A 



1 



4A, 



(26) 



To gain a sense of the magnitude of the mass matrix 
elements, we assume the following values for the observed 
mass-squared differences: 

A L » 1 eV 2 , 

A A ps 3 x 10~ 3 eV 2 , 

A s » 10~ 5 eV 2 , (27) 

and so the ratios of mass-squared differences are all the 
same, namely 



A A As 

^i«_l«3x 10" 
A L A a 



(28) 



It is interesting to note that, for the above value of Al, 
this is also the value of the weak mixing angle between 
upper and lower doublets needed to fit the LSND data 



sin 2 289 ps 3 x 10" 



(29) 



The large parameter in the mass matrix, tuq, is close 
to 0.5 = eV: 



2ron 



1.001 eV 



(30) 



and the small parameters, e and 2rrid, are much smaller 
and roughly equal to one another, 



e w 1.5 x 10~ 3 eV , 2m d ps 1.5 x 10" 3 eV 



(31) 



Thus the upper doublet, corresponding to v T and v^, 
has a central value of 1.001 eV and a splitting of 1.5 x 
10~ 3 eV, while the lower doublet, corresponding to v e 
and u s , has an almost zero central value, 3 x 10 -3 eV, 
with the same splitting as the upper one. 

We have not considered the Majorana versus Dirac na- 
ture of the four neutrinos and the constraints from no- 
neutrino double beta decay If the three active ones 
are all Majorana particles, then the sum of their masses 
times CP phase must not exceed the current bound of 
0.2 — 0.6 eV jL4|. In the above example, this is most easily 
achieved by giving the members of the upper doublet op- 
posite CP phases, which make them "Pseudo-Dirac" neu- 
trinos because of the small mass difference 2ra<j. What- 
ever phase is assigned to the active member of the lower 
doublet, the sum of masses times phase will not exceed 
6 x 10~ 3 eV, well within the experimental limit [ fL5| ]. 

We may now ask whether the mass matrix M4 can be 
derived from a symmetry principle. As we have noted 
above, the case of maximal mixing among the two mem- 
bers of a doublet corresponds to the permutation symme- 
try S2 between them. Likewise the general structure of 
M4 involves the permutation symmetry S2 between the 
two doublets. It is not difficult to show that the most 
general 4x4 matrix which is invariant under S% x S2 
is given by: 



Hi 



X Y 
Y X 



(32) 



The (2 x 2) submatrices X , Y are both of the same S2 
symmetric form as M2 above. 

Comparing M4 with H4, we see that it is of exactly the 
same form except that the off-diagonal submatrix K is a 
multiple of the unit (2 x 2) matrix whereas Y can have 
an off-diagonal matrix element. Physically, the absence 
of an off-diagonal matrix element in K means that the 
splitting between the members of the upper doublet is 
exactly the same as that between the members of the 
lower doublet — an "equal splitting" rule. 

In conclusion, we have constructed a mass matrix 
which can simultaneously accommodate all three indi- 
cations for neutrino oscillations. Its particular structure 
as a direct product of (2 x 2) matrices can be derived from 
an underlying 52 x S2 symmetry plus an equal splitting 
rule. It may be interesting to speculate that this symme- 
try might in turn be a subgroup of a larger permutation 
symmetry, for example S4 , and that the larger symmetry 
can be used to distinguish between the active and sterile 
neutrinos. For example, the three active neutrinos could 
belong to a triplet with respect to an S3 subgroup of the 
larger group, while the sterile neutrino is a singlet. 

We recognize that large mixing between a sterile neu- 
trino and the electron-neutrino in the Solar Neutrino 
Problem can disturb Big Bang Nucleosynthesis, |T(| and 
we have no ready solution for this problem. Whether 
Big Bang Nucleosynthesis can accommodate 3 or 4 light 
neutrino degrees of freedom will depend crucially on the 
amount of primordial deuterium in the universe; at the 
moment this is not well determined. jl7) We do, how- 
ever, regard the existence, or non-existence, of a sterile 
neutrino to be an experimental question which will even- 
tually be settled by the observation of the neutral-current 
interactions of solar neutrinos, as in the SNO experiment. 



We are indebted to Hamish Robertson for asking a 
question which sparked this investigation. 
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